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C"| ' Abstract. The conversion of an ordinary wave to an extraordinary wave in a 2D 

inhomogeneous slab model of the plasma confined by a sheared magnetic field is studied 
analytically. 
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The linear conversion of waves in nonuniform media is a common phenomenon in physics 



and has been studied in fields as diverse as RF heating of fusion plasmas, ionospheric 
physics and so on. The particular case of the linear mode conversion, namely, the 
\ conversion of an ordinary (O) wave to an extraordinary (X) wave at electron cyclotron 
frequencies being a primary consideration for an auxiliary electron heating and current 
drive in fusion devices was intensively studied for plasmas with a cold plasma dielectric 
tensor. A quantitatively accurate description of this problem was paid much attention to 
a couple of decades ago within both the WKB approximation and the full wave analysis 
for the ID inhomogeneous slab model of the plasma with zero shear pQ- [3]. Having been 
renewed by Wendelstein 7-AS (W7-AS) considerable success in realization of the mode 
conversion scheme [5] the problem was examined in a number of papers [6], [7j, [8] within 
the 2D plasma model and in [9] within the 3D plasma model by analysis of a set of the 
reduced wave equations valid in a vicinity of the O-X mode conversion layer. In all these 
papers the curvature of the magnetic flux surfaces was neglected under an assumption 
of high localization of the conversion region and magnetic field was considered to be 
shearless with the straight magnetic field lines. No effort was mounted to describe the 
O-X mode conversion in more realistic case of the 2D inhomogeneous plasma in the 
magnetic field with a shear except possibly the paper [10] where the effect of the shear 
on the efficiency of the O-X mode conversion was studied within the ID inhomogeneous 
slab model of the plasma. Here we are concerned with the analysis of the O-X mode 
conversion problem for the 2D inhomogeneous plasma with the sheared magnetic field. 
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2. Introduction 

The linear conversion of waves in nonuniform media is a common phenomenon in physics 
and has been studied in fields as diverse as RF heating of fusion plasmas, ionospheric 
physics and so on. The particular case of the linear mode conversion, namely, the 
conversion of an ordinary (O) wave to an extraordinary (X) wave at electron cyclotron 
frequencies being a primary consideration for an auxiliary electron heating and current 
drive in fusion devices was intensively studied for plasmas with a cold plasma dielectric 
tensor. A quantitatively accurate description of this problem was paid much attention 
to a couple of decades ago within both the WKB approximation and the full wave 
analysis for a slab model of the ID inhomogeneous plasma with zero magnetic field 
shear []]- [3]. Having been renewed by Wendelstein 7-AS (W7-AS) considerable success 
in experimental realization of the O-X mode conversion scheme [5] the problem was 
examined in a number of papers [6]- [?] within the 2D plasma model and in [9] within 
the 3D plasma model by analysis of a set of the reduced wave equations valid in a 
vicinity of the O-X mode conversion layer. In all these papers the curvature of the 
magnetic flux surfaces was neglected under an assumption of high localization of the 
conversion region and magnetic field was considered to be shearless with the straight 
magnetic field lines. No effort was mounted to describe the O-X mode conversion in 
more realistic case of the 2D inhomogeneous plasma in a sheared magnetic field except 
possibly the paper [10J where the effect of the shear on the efficiency of the O-X mode 
conversion was studied within the ID inhomogeneous slab model of plasma. In this 
work we are concerned with the analysis of the O-X mode conversion problem for the 
2D inhomogeneous plasma with the sheared magnetic field. 

3. Basic equations 

In this paper we consider plasmas with 2D inhomogeneity that assumes the surfaces 
n = const and B = const being not approximately parallel and the external sheared 
magnetic field. We treat a problem in the simplest slab geometry accounting for plasma 
axisymmetry, expecting that the results being obtained within this simple model hold 
in toroidal plasmas. 

Let us start with introduction of the Cartesian co-ordinate system (x, y, z) with 
its origin located at the O-mode cut-off surface and coinciding with the center of an 
incident beam of the monochromatic ordinary waves. The co-ordinates x, y and z being 
scaled in c/u, where a; is a frequency of the wave and c stands for the speed of light, 
imitate the flux function, the poloidal angle and the toroidal angle, respectively. We 
define the magnetic field as 

B = Bq (sin 6e y + cos 6e z ) , 

where 9 = 8q + 68, 9q is an angle between eii = B/5o and e z at x = and 56 = 56(x) 
(69 = at x = 0) originated due to the shear of the magnetic field. In this study we 
assume for simplicity that 8q = and en \\e z at x — 0. The parameters of inhomogeneity 
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are assumed to be small |Vln (n)| <C 1 and |Vln (B) \ -C 1 as in the previous studies 0- 
[9]. In such a case the WKB approximation is applicable for description of the interacting 
waves except a narrow localized conversion layer in a vicinity of the O-mode cut-off 
surface where the full-wave analysis connecting the incoming and outgoing WKB waves 
is needed. The WKB approach can give us not only the boundary conditions but 
some hints on how to reduce the full-wave equations in the conversion layer as well. 
The dispersion curves of the O and X modes are adequately described by the cold 
plasma model. These curves are coupled if the local wave-vector n with components 
n x = n ■ e X} n v = n • (e\\ x e x ), n\\ = n • ey has at the O-mode cut-off surface the only 

parallel component equal to n res = (1 + u /uo ce ) where u ce is the electron cyclotron 
frequency at x — 0, y — 0. This means the conversion of the O wave into the X wave 
occurs if in a close vicinity of the O-mode cutoff surface n\\ ~ n res + <5n|| where 5n\\ <C 1 
and n x -C 1, n v <C 1. We expand the components of the cold dielectric tensor and 
the angle 9 near the origin of the co-ordinate system into the Taylor serious up to the 
first-order terms: 

e + — n res ~ — |Ve+| (xsimp — ycosy?) , e\\ ~ — |Vq||x, 9 ~ Gx, 

where ip is an angle between two surfaces e + — r? Tes = and en = 0, G' 1 is a length 
describing the magnetic field shear. The wave vector in the conversion layer can be 
interpreted as the differential operator 

n x = —id x , n r] = —id y + Gx, m\ — n res = id z , d r = d/dr 

acting on an appropriate component of the electric field. Using the standard expansion 
procedure over the small parameters and keeping the first order terms one can obtain 
the following set of equations for the electric field components 

n res D + E z — i (|Ve+| (a; cosy? — ysmip) — i2n res d z ) E + = 

i2\ Vey \xE z - n res D_E + = (1) 

E_ = 0, 

where E± — E x ± iE y , D± = d±^f n res Gx, d± = d x ± id y . Variation of the electric field 
components in (1) is neglected as it results in the effects being higher order of vanishing. 
For any detail of (1) derivation we can refer readers to [?]. This set of equations is 
different from that studied in [?], which describes the 2D O-X mode conversion but 
misses the term n res Gx originated due to the magnetic field shear. 

This is our basic system which as it stands lost symmetry inherent in the case of 
the magnetic field with zero shear G = 0. It is at first disturbing us by virtue of the fact 
that no separation of variables and no existence of eigenvalues are permitted in the case 
of the sheared magnetic field. Fortunately, this initial impression is misleading as will 
be shown in the forthcoming sections. With the use of the procedure proposed in [?] 
we will obtain the second order partial differential equation instead of the system (1). 
Then, applying an integral representation of Laplace integral type as in [?] of unknown 
function, we will separate variables and obtain eigenvalues and eigenf unctions of this 
partial differential equation. 
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4. Analysis of the system (1) 

Following a procedure having been proposed in [?] opens a way of overcoming underlying 
mathematical difficulties and permits us to separate variables in (1). As we assume an 
axisymmetric plasma, we can seek a solution of (1) in a form ~ exp (in z z), where 
n z = const. Introducing new notations 

2|V5|||/|V £+ |) 1/4 ^, (2) 



V5 + |/(2|Vq||)) 1/4 £ + , 

g = Gnlj(2\Ve ll \\Ve + \) 1/2 , 
y - 2n res n z / smip/\Ve + \ -> y, 

x,y (2|V£||||V£+|) 1/4 /^rL 2 -*x,y 

we read (1) as 

iD + a + (x cos ip — y sin ip) b = (3) 
iDJb + x ■ a = 

where D± = d± =F g% and d± = d x ±id y . Let us use a generalized Laplace transforms of 
a(x, y) — ► a(x, q) = L [a(x, y)\ and b(x, y) — ► 6(a;, q) — L [b(x, y)\ where 

^ [•] = exp \-%-r— Q 2 I / rfyexp (-iqy -i^^y 2 ) [■] (4) 
\ smip J Jc \ 4g J 

As we will see later, a contour C can be chosen along the real axis of the complex 

variable q. We introduce k = —2g/ sin<^ • q and seek a(x, k), b(x, k) in the form 

a(x, k) — d-f — ((g — i/3 cos cp)x + if3 sin <pk + 2/3gdk) ■ f, (5) 
b(x,k) = P-d + f + (g/3 + i)x-f, 

where d± = d x ± idk, — (£ — exp ^ = (cos (/? — g 2 — i sin <^) 1//2 = |£| exp (—iip), 
\C,\ = (1 — 2(7 2 cos if + g 4 ) 1 ^, 2ip = arctan [sin (p/ (cos ip — g 2 )]. Using (4) and (5) we 
derived a partial differential equation of second order for a function / instead of the 
system (3) of two partial differential equations of first order for a and b 

d lf + d lf + ((cos^-# 2 )x 2 -mn<p-xk + i\Z\exp(-ii/;)) f = (6) 

This equation deserves few comments. For zero shear, g = 0, (6) coincides with the 
equation analyzed in [?]. For the sheared magnetic field we estimate the parameter g 
as g ~ L n j L Q 1, where L n is a scale on which the the density profile at the O-mode 
cut-off surface changes and L is a minor radius. Henceforth, we keep in mind that the 
coefficient cos <p — g 2 at x 2 in (6) is positive. At tp = we return to the equation treated 
within the ID O-X mode conversion problem in the sheared magnetic field. First who 
analyzed this case by the WKB analysis were Cairns and Lashmore-Davies pi)] . To 
separate variables in (6) we introduce 

u = \f\£\ (cos ip • x — sin ip ■ k) 

|£| (sin?/> • x + cosip ■ k) 
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that yields 

dlf + d 2 J + (cos^V - sin^V + iexp (-#)) / = (7) 
A particular solution of (7) is 

oo 

f(u, v) = J2 c P Diu P h ( u ) ( v )> ( 8 ) 

p=0 



where u p = 7r|tan^| (p+(l - sign(ip)) /2), D iVp/lT (u) = D iUp/lT [y/2 cosipexp (m/A)u 

is the parabolic cylinder function [IT] . p (t>) = P (^/sin \ib\vj are the Hermitian 
polynomials 

(j> p {v) = (2^P ] Y 1/2 exp(-v 2 /2)H p (v) 

possessing a property / (bp (v) <bk {v) dv = 5 p k- Substituting (8) in (5) we obtain 

a(x, k)=R (exp(#) |£| + ig) - I (exp(#) \t\ - w) (9) 
b(x,k) = R + I, 



where 



1 = V If I c P D iv P fr ( u ) ( d v + sin if) ■ v) <p p (v) , (10) 

p=0 
oo 

R = yl^lY. ( v ) {~ id u + cosip ■ u) D iUp/n (u) . (11) 

p=0 



Simple but tedious algebra yields 



/ = yj\ sin(^)| £ C P A>/. («) P («) , (12) 

p=0 



R = \Ji cos(^) 5^ C p A 7p /7r-i («) A p (v) , (13) 

p=0 

A p = 4> p+ i (v) at ip > 0, A p = —2pcf)p_i (v) at tp < 0, 

where 

7 P = ir\t8Jiif)\(p + ^(l + sign((p))\. (14) 

The term J corresponds to a wave propagating in the positive direction along u, 
namely to the incident wave and converted wave. The term R being proportional to 
Di 7 / n -i(u) corresponds to the reflected wave propagating in the negative w-direction. 
The coefficients C p in (12) can be chosen in such a way to fit the incident WKB ordinary 
wave in the WKB region. To define C p we perform the transform (4) of the incident 
beam A (x, y) — > A (x, q) = L [A (x, y)]. Using the definition of k, u and v we express 
A(x,q) in terms of these variables A(x,q) — ► A(u,v). We assume A(u,v) at u b — > — oo 
in the form 

/ COS lb 'Y / I \ \ ~ 

A(u, v) = exp lin res z — i — - — u\ + i— In f y 2 cos ipUbj \A(v), 
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As far as (12) is an expansion in orthogonal functions (j) p , we can define C p as follows 

/oo 
dvA{v)(j) p {v) (15) 
-oo 

What distinguishes (15) from the coefficients obtained earlier in [?] for the magnetic 
field with zero shear is the integrand A(v ) which is wider function of v due the transform 
(4) than in the case g = 0. As a consequence of this the eigen-mode spectrum C p is 
changed considerably compare to the case of the magnetic field with zero shear. 



5. Conversion coefficients 



Since I corresponds to the incident wave at u < and to the converted wave at u > 
we can derive the conversion coefficient for the O-X mode conversion using asymptotic 
expression of Di yp / n (u) [H] entering (12). Let us introduce the conversion coefficient 
(T p ) corresponding to a separate p term in the sum (12) over the Hermitian polynomials 

T p = exp (— 2n(p + fi) tan |-0|),/x = 1 at 93 > 0, fi = at (p < (16) 

2 \ 1/2 

, siny? 2 / (cosy? — g 2 ) 

tan \ip I 



! 1 + siny? 2 / (cos<^ — g 2 ) 2 + 1 

As in the case of the magnetic field with zero shear, the conversion coefficient for the fixed 
p mode (16) has asymmetry, considered first time in [?], with respect to the sign of the 
angle p between two surfaces e\\ = and e + — n 2 es = (// = 1 at ip > 0, // = at p < 0). 
For the fixed n z the conversion coefficient T for the O-X mode conversion may be 
expressed as 

T= E,T\C/ 

For the magnetic field with zero shear g = the transform (4) yields A (x, y) = 
L [A (x, y)\ and the conversion coefficient of the fixed mode p is performed as in [?] 

T p = exp (-27rtan (\<p\/2) (p + //)). 

At p = with the use of Mehler's formula [TT] we return to the conversion coefficient 



T 



Jdq 


A(q) 


2 exp (-2vrg 2 / (1 " 9 2 f 2 ) 


fdq 




2 



A(q)= J dy- A (y) exp (-iqy) 



derived firstly by Cairns and Lashmore-Davies in [10J by the WKB analysis. We remind 
that q and g are scaled according (2) and n z = is assumed. Analyzing (4), (9), (12) and 
(14) we can summarize that the magnetic field shear degrades the efficiency of the O-X 
mode conversion. An important feature is increasing of the degradation with increasing 
of the mode number. Though the parameter g for the typical experimental conditions 
is small g ~ L n /L C 1 we can expect remarkable effect of the magnetic field shear on 
the efficiency of the O-X mode conversion for the modes with p> 1. 
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6. Conclusion 

In the paper new insight on the 2D O-X mode conversion problem has been done by 
considering plasmas confined by the sheared magnetic field. Ignoring the curvature of 
the magnetic flux surfaces under an assumption of high localization of the conversion 
region and assuming the straight magnetic field lines we have studied the problem in the 
simplest slab geometry accounting for plasma axisymmetry. Nevertheless, the results 
being obtained is expected to hold in more realistic toroidal plasmas. The resulting 
expressions for the electric field components (9), (12) and (13) have been derived that 
permits finding the conversion coefficient (16) explicitly. As it has been demonstrated 
the magnetic field shear degrades efficiency of the O-X mode conversion. For usual 
experimental setup this degradation appears to be remarkable. Finally, extending the 
model to the case when the magnetic field direction do not parallel to the toroidal 
direction at the O-mode cut-off surface which is more relevant to tokamak physics needs 
to be done. This issue will be studied in future paper. 
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